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Abstract. Nuclear structure data are of crucial importance in order to address impor- 
tant astrophysical problems such as the origin of chemical elements, the inner working 
of our Sun, and the evolution of stars. We demonstrate this by investigating the ground 
state structure of 8 B and 7 Be nuclei within the Skyrme Hartree-Fock framework and 
' by calculating the overlap integral of 8 B and 7 Be wave functions. The latter is used to 

£N) calculate the astrophysical S factor (S17) for the solar fusion reaction 7 Be(p,7) 8 B. 
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1. Introduction 
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Nuclear astrophysics research addresses some of the most fundamental questions 
in nature e.g., the origin of the elements that make up our bodies and our world, 
and formation and evolution of the sun, the stars and the galaxies. There is an 
intimate connection between nuclear physics inputs and studies of these fascinating 
astrophysical phenomena [1]. A diverse set of nuclear data is required to model 
the composition changes and energy release in astrophysical environments ranging 
from Big Bang to inner working of our own Sun to exploding stars. Theoretical 
studies as well as measurements of microscopic nuclear physics processes provide 
the foundation for models to understand various astrophysics phenomena. These 
models are now challenged by incredibly detailed observations from ground based 
and space based (e.g., CHANDRA X-ray Observatory, Hubble Space Telescope) 
observational devices that give us an unprecedented view of the Cosmos 

We shall demonstrate the role of the nuclear structure physics in nuclear astro- 
physics by considering the case of the Solar fusion reaction 7 Be(p,7) 8 B, which is 
of great importance to the solar neutrino issue and to other related astrophysical 
studies. 8 B is the source of the high energy neutrinos from the Sun that are detected 
in the SNO, Kamiokande and Homestake experiments [2, 3, 4]. It is, therefore, a 
crucial requirement to determine as accurately as possible the cross section of this 
reaction at relative energies corresponding to solar temperatures (about 20 keV). 
In this energy region, the cross-section a pl (E cm ) [which is usually expressed in 
terms of the astrophysical Sn(E cm ) factor] of the 7 Be(p,7) 8 B reaction is directly 
proportional to the high energy solar neutrino flux. A better knowledge of S17 is, 
therefore, important to improve the precision of the theoretical prediction of 8 B 
neutrino flux from present and future solar neutrino experiments. 
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617(0)] is determined cither by direct measurements [5, 6] or by indirect methods 
such as Coulomb dissociation [7] and transfer reactions [8, 9]. Efforts have also been 
made to calculate the cross section of this reaction within the framework of the 
shell model and the cluster model [10, 11]. The key point of these calculations is 
the determination of the wave functions of 8 B states within the given structure 
theory. 

We have investigated the structure of 8 B in the framework of the Skyrme Hartree- 
Fock (SkHF) model which has been used successfully to describe the ground-state 
properties of both stable [12, 13] as well as exotic nuclei [14]. The SkHF method with 
density-dependent pairing correlation and SLy4 interaction parameters has been 
successful in reproducing the binding energies and rms radii [14] in the light neutron 
halo nuclei 6 He, 8 He, n Li and 14 Be. We solve spherically symmetric Hartree-Fock 
(HF) equations with SLy4 [15] Skyrme interaction which has been constructed 
by fitting to the experimental data on radii and binding energies of symmetric 
and neutron-rich nuclei. Pairing correlations among nucleons have been treated 
within the BCS pairing method. We have, however, renormalized the parameter 
of the spin-orbit term of the SLy4 interaction so as to reproduce the experimental 
binding energy of the last proton in the 8 B nucleus [16]. A check on our interaction 
parameters was made by calculating binding energies and rms radii of 7 Bc, 7 B, 
8 Li and 9 C nuclei with the same set where a good agreement is obtained with 
corresponding experimental data. The overlap integral of the HF wave functions 
for 7 Be and 8 B ground states has been used to calculate the astrophysical Sir 
factor. 

2. Results and discussions 

2.1. Structure calculations 

The values of various parameters of the SLyA Skyrme effective interaction as used 
in our calculations are given in [16]. The rms radii for matter (r m ), neutron (r„) and 
proton (r p ) distributions are presented in Table I for five light nuclei. Also shown 
in this table are the matter, neutron, and proton radii (under the column "Expt") 
extracted by methods in which measured reaction (or interaction) cross sections are 
fitted by theoretical models having them as input parameters. The quantities listed 
under "Theory" column are the results of our calculations. Here r m is obtained 
by summing the average of proton and neutron radii in every orbit weighted with 
occupation probabilities. We see that for all the isotopes the calculated r m is in 
good agreement with the corresponding values listed under the "Expt" column. 

2.2. Valence proton radius in 8 B and Astrophysical Sir factor 

We define a overlap function of the bound state wave functions of two nuclei B 
and A, where B = A + p (p represents a proton) as 




(1) 
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Table 1: Rms mass (r m ), proton (r p ) and neutron (r n ) radii for various nuclei. 
Theoretical results obtained with modified SLy4 Skyrme force. Under the column 
"Expt", are the values of corresponding radii extracted by fitting the reaction or 
interaction cross sections by different theoretical methods as discussed in the text. 
We have defined n —< rf > 1 / 2 



Nucl. 




rms 


radii (fm) 












"Expt" 






Theory 








r P 


r n 




r P 


r n 


7 Be 


2.33 ±0.02 






2.49 


2.63 


2.29 


7 B 








2.86 


3.18 


1.84 


*Li 


2.37 ±0.02 


2.26 ±0.02 


2.44 ±0.02 


2.54 


2.29 


2.67 


S B 


2.55 ±0.08 


2.76 ±0.08 


2.16 ±0.08 


2.57 


2.73 


2.27 




2.43 ±0.03 


2.49 ±0.03 


2.33 ±0.03 








9 C 


2.42 ±0.03 






2.59 


2.77 


2.20 



where Ia and Ib are the total spins of nuclei A and B, respectively, and r is the 
position of the proton with respect to the cm. of nucleus A. a p is the spin variable 
of the proton, and £ stands for the remaining set of internal variables which also 
include isospins. In this expression nuclear wave functions \I> J 4 and tyg are supposed 
to be properly translational invariant. The Hartree-Fock wave functions calculated 
in the previous section may not be so despite the fact that a cm. correction factor 
has been incorporated in the energy functional. Corrections for the spurious cm. 
motion should, therefore, be applied to the HF wave functions before using them 
in Eq. (1). However, effects of such corrections on the one-particle overlap function 
calculated within the shell model have been found [17, 18] to be of the order of only 
about 2-5 %. It is quite likely that the situation will be no different for the HF case. 

We, now, present results for the astrophysical Si 7 factor calculated using HF 
overlap functions. In the region outside the core where range of the nuclear inter- 
action becomes negligible, the radial overlap wave function of the bound state can 
be written as 

RMi(r) - cyW v . l+1/2 (2kr)/r, (2) 

where W is the Whittaker function, k the wave number corresponding to the single 
proton separation energy and 77 the Sommerfield parameter for the bound state. 
In Eq. (2), c7j is the asymptotic normalization constant, required to normalize the 
radial overlap wave function to the Whittaker function in the asymptotic region. 
The S17 factor is related to the proton capture cross section as 



S 17 (E) = a{E)E e^ E) l 



(3) 
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At the zero energy the Si 7 factor depends only on Cij and one can write [10, 11] 

S ir (0) = k^^j, (4) 

3 

In Ref. [10], k (= 37.8) has been obtained by using a microscopic cluster model 
for the scattering states, while a value of 36.5 has been reported for this quantity 
in Rcf. [11] using a hard sphere scattering state model. However, once the relevant 
integration distances are sufficiently enhanced in [11] the value of k there comes 
out to be 37.2, which is in good agreement with that of the microscopic model. 

In our calculations, we have used our HF overlap function directly in the cal- 
culation of the amplitudes for the direct capture reaction. The amplitude for the 
radiative capture reaction can be written asfe + c^a + 7 which can be written as 

M fi = {M&MidWM&MtcWSHr)), (5) 

where <fr, £ and r are the wave function, and the internal coordinate of the bound 
particle and the relative coordinate between b and c, respectively. 0(r) is the elec- 
tromagnetic operator, ip^ (r) is the distorted wave in the initial b + c channel. The 
overlap of initial and final bound state wave functions is defined by 

7 b a c (r) = r)|&(&)&(&)>, 

= ^ < I b M b jm\I a M a >< Im - \isii\jm > i l B% tj {r)Y tm - ,,,(?) 

x x»„, (6) 

It may be noted that at low relative energies the distortion effects are goverened 
solely by Coulomb interactions. Substituting Eq. (6) into Eq. (5), expanding ^ + ^(r) 
in terms of the partial waves and carrying out the angular momentum algebra, the 
amplitude Mfi can be written in terms of the wave function R^- (r) and the radial 
coulomb wave functions for the initial channel. Therefore, the overlap functions 
Rb£j(r) calculated in the HF theory can be used to calculate the capture cross 
sections at low relative energies without any uncertainty of other input quantity. 
It may further be noted that using Eq. (2), we can also calculate the asymptotic 
normalization constants ctj. 

In Table II, we show the results for the asymptotic normalization constants for 
the P1/2 and p 3 / 2 proton orbits and the the astrophysical 5*- factor obtained by this 
method. Our c 3 / 2 agrees almost perfectly with the value of this quantity determined 
experimentaly in Ref. [19]. However, our C1/2 is larger than that reported by these 
authors by a factor of about 1.5. They have determined these quantities (by invoking 
mirror symmetry) from the ANCs of the 8 Li — > 7 Li +n system extracted from the 
measurement of the neutron transfer reaction 13 C( 7 Li, 8 Li) 12 C. This difference in 
the calculated and "measured" values of C\i<i is significant as a larger C\ii leads to 
a bigger 5 17 as compared to that reported in Ref. [19]. There may be a need to 
relook in the theoretical analysis of the data reported in Ref. [19]. In the DWBA 
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calculations of the 13 C( 7 Li, 8 Li) 12 C reaction, same set of potentials have been used 
for 7 Li + C and 8 Li + C channels. However, 8 Li is an unstable system; thus optical 
potential in the final channel could be different from that of the incoming ( 7 Li) 
channel. It is also not very clear from Ref. [19] if the spin-orbit terms have been 
used in the optical potentials and bound state potentials. This could affect the pi/ 2 
and p 3 /2 components differently. More work is clearly needed to clarify this issue. 



Table 2: SkHF results for asymptotic normalization coefficients (<%), and astrophys- 
ical S— factors 5^(0) and 5^(0) and Si 7 (0). Sf 7 corresponds to results obtained 
by using the HF overlap function directly to a capture code while Si 7 and S^ 7 to 
those obtained by using Eq. (4) with k = 36.5 and 37.8, respectively. 



Observable Our Result 



C13/2 0.64 

en/2 0.34 

Sf 7 (0)(eVb) 22.0 

Sf 7 (0)(eV6) 19.5 

S? r (0)(eVb) 20.2 



The recent 7 Be(p, 7) measurements yield values of Si7(0) which are clustered 
around 18.5 eV b [5] and 22.0 eV b [6]. Our 5i 7 (0) (22.0 eV barn) is closer to two 
latest direct capture measurement results which claim good accuracy. However, this 
is larger than the most recent result reported by the Texas A & M group from the 
ANC method [19]. This is also slightly larger than the values obtained from the 
Coulomb dissociation (CD) method [20], even thought it has been argued [21] that 
there is no significant difference between the CD and diect capture values of £17(0). 



3. Summary and conclusions 

In summary, in this paper we studied the structure of 8 B, and 7 Be nuclei within 
the Skyrme Hartree-Fock (SkHF) framework. We calculated binding energies, var- 
ious densities distribution and rms radii for these nuclei. Using the same set of the 
force parameters, we obtain good agreements with experimental values of binding 
energies and rms matter radii for all these nuclei. We have calculated the overlap 
function < 7 Be| 8 B> from the SkHF wave functions which has been employed to 
extract the asymptotic normalization coefficients for the 8 B ^ 7 Be +p system. We 
obtain an astrophysical S-factor of 22.0 eV b which lies within the adopted limits 
(19.1^i'o e V b) of the near zero energy astrophysical S-factor. Our work clearly 
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shows that proper nuclear structure input is vital in order to understand some very 
important issues in nuclear astrophysics. 
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